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Flavour singlet physics in lattice QCD with background fields 
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We show that hadronic matrix elements can be extracted from lattice simulations with background 
fields that arise from operator exponentiation. Importantly, flavour-singlet matrix elements can be 
evaluated without requiring the computation of disconnected diagrams, thus facilitating a calculation 
of the quark contribution to the spin of the proton and the singlet axial coupling, g° A . In the two 
nucleon sector, a background field approach will allow calculation of the magnetic and quadrupole 
moments of the deuteron and an investigation of the EMC effect directly from lattice QCD. Matrix 
elements between states of differing momenta are also analysed in the presence of background fields. 



The composition of the spin of the proton is a question of long-standing interest since experimentally it is known 
that only about 20% of the total spin of the proton is carried by the quark helicity 1] . Contributions from quark 
orbital angular momentum and from gluon angular momentum are obviously important, but are currently unknown. 
Recently Ji |2j has shown that a well-defined decomposition of the proton spin exists, 



2 - Jq + J 9 - g 



J q + Jg = -E, + L q + Jg , (1) 



where each term is gauge invariant (though scheme and scale dependent) and given by proton matrix elements of the 
operators 

J q = [d 3 rrxT q (2) 
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3 g = J d 3 rr x (E x B) , (3) 

where colour indices are suppressed, = T|° and T£ v = ipj^iD^ip is the traceless, symmetric part of the quark 

energy-momentum tensor (D — |[£) — D]). The quark helicity contribution, S g , has been measured in polarised 
deep-inelastic scattering, but the separation of the quark orbital angular momentum and the total gluon angular 
momentum has not yet been determined from experiment though it can be extracted from deeply-virtual Compton 
scattering. 

A reliable prediction of the various contributions to Eq. Q from QCD would obviously be an important achievement. 
However, in QCD this question is non-perturbative since it involves hadronic-scale physics, and to address the issue 
from first principles, one must use lattice QCD. Indeed,proton matrix elements of both the helicity operator, ip^^ip, 
and the energy momentum tensor, T£ v , have been computed and some estimates of the partitioning 

of the proton spin have been made. These calculations have been performed at unphysically large quark masses 
and on modest volumes and such parametric limitations make the connection of these results to the physical world 
non-trivial, though progress continues to be made in this regard 0, IToL [Til IT2I IT^ | . Such issues must be confronted for 
all hadronic quantities extracted from lattice simulations. However, the total quark helicity (singlet axial coupling) 
and angular momentum content of the proton are flavour singlet quantities {g\ ~ S g ~ S M + Ed and J q ~ J u + J4) 
and lattice calculations of them are further plagued by so-called quark-line disconnected contributions in which the 
relevant operator is inserted on a quark line connected to the proton source only by gluons as in Fig.^b). These terms 
are notoriously difficult to compute (see Ref. |l4j for a review) and result in relative errors an order of magnitude 
larger than the corresponding connected contributions. 

Numerous other quantities such as the strangeness magnetic moment and the pion-nucleon sigma term suffer from 
the same difficulties. In the two nucleon sector, notable examples are the deuteron magnetic and quadrupole moments, 
and the deuteron structure function, F${x), (accessible in lattice QCD only through two-particle matrix elements of 
twist- two operators ^|) that gives the simplest manifestation of the EMC effect [lr^ . Even were it feasible to 
calculate five point correlators on present day computers, these quantities would again require orders of magnitude 
more computational effort than their flavour non-singlet analogues. 

In this article, we highlight an alternative method to calculate operator matrix elements that eliminates the issue of 
disconnected contributions at the expense of requiring the generation of additional ensembles of gauge configurations. 
This essentially involves computing two-point correlators in the presence of generalised background fields which we 
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FIG. 1: Connected (a) and quark-line disconnected (b) contributions to a baryon matrix element. The squares correspond to 
the baryon source and sink and the dark circle to the operator insertion. 



show arise from exponentiating the operator whose matrix element we wish to calculate into the QCD action (see 
Ref. |17| for an early example). Since disconnected contributions are included automatically, this approach may 
provide a means of extracting reliable results for singlet matrix elements. 

Electroweak background field methods have a long history in lattice QCD. In the early 1980s, Martinelli et al. 
|l8j | and Bernard et al. |l9j | performed the first calculations of the magnetic moments of the proton, neutron and 
A-resonance by calculating the baryon masses on a lattice immersed in a constant magnetic field. This approach has 
since been extended to calculate electric and magnetic polarisabilities of various hadrons |2Ctl2ll l22. 23]. Fucito et al. 
[l7| also showed that a background axial field could be used to determine the isovector axial coupling of the nucleon, 
gA- Recently it has been demonstrated [24| that background electroweak fields can be used to probe the electroweak 
properties of two-nucleon systems. Measurement of two particle energy levels at finite volume in magnetic and weak 
background fields is sufficient to determine the weak- and magnetic- moments of the deuteron, and the threshold 
cross-sections for radiative capture (np — > d'j) and weak-disintegration of the deuteron (e.g. v x d — > v x np [x = e, [i, r] 
and v e d — * e + nn). 

In this work, we extend these analyses to consider more general background fields that correspond to exponentiation 
of a larger class of operators. We additionally consider the use of background fields to extract operator matrix elements 
between states of differing momenta. Such matrix elements determine the electromagnetic form-factors and moments 
of generalised parton distributions. Whilst much of the following analysis will be equally applicable in the pure-gauge 
and mesonic sectors, we restrict our discussion to baryons. 

One baryon: In the single baryon sector, one is interested in calculating matrix elements of the form 

(B'{p', s')\0\B(p, «)) ~ u B ,(p', S ')Mb>b[0}u b (p, s) , (4) 

where q — p' — p is the relative momentum of the baryons B and B' (s and s' label their spin) and ub{p, s) is a 
baryon spinor. The Dirac structure Mb'b[0] is parameterised by a set of scalar form factors, Fi(q 2 ), depending on 
the particular operator, 0, and the external states (for the case of the energy-momentum tensor, we shall be more 
specific below). Typically, such form factors are calculated by computing ratios of two and three point correlators. 
We define 

G BB (p,t;T) = Y,e- ip - x T Pa (0\xi B H^t)x { f } (0,0)\0) , (5) 

X 

and 

G B , B(p,p',t,T;T') = J2 r^e^ x e l P'- z (0|xi B ' ) (x,i)O(y I r)x? ) (z,0)|0), (6) 

where T and V are Dirac projectors (the exact form of which depends on the operator under consideration), and 
Xq S ' ) (x^) i s an interpolating field for the baryon sink (source); for the proton, a common choice is Xa\x) = 
e ijkUa(x)u : 'p(x)(C^5) pyd^x) . It is easy to show (by inserting complete sets of states, see e.g. Ref 0) that the ratio 

G BB (p / , t; T)G bb (p', t sink ; T)G BB (p, t sink - r; T) 1 1/2 
G BB {p,r;r)G BB (p,t sin y,r)G BB (p',t siBk -T;T)\ ' U 

will exhibit a plateau over a range of time slices a <^ t <ti £ s i n k where the sink time, tsink, is held fixed. For time 
slices within this plateau and appropriate choices of the Dirac projectors, this ratio is proportional to the form-factors 



H(p , p, t sink , T,l,L ) - — — — 

(-f-B-BlP j tsink, t J 
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under consideration. This approach and its variants (see e.g. Refs. H^,[2^|) are referred to as the operator insertion 
method. For flavour singlet operators, calculation of Gbob requires evaluation of both types of diagram in Fig. ^ 

For a local operator, O(z), at most bilinear in quark fields, it is possible to calculate the same matrix element by 
measuring a two-point correlator on an ensemble of field configurations generated with a modified Boltzmann weight. 
To see this, we add a term 

d d zO(z) -n(z), (8) 

to the usual QCD action, Sqcb, with a fixed external source field fl(z). If O(z), which depends on the quark (ifj) 
and gluon (A^) fields and their derivatives, carries Lorentz or flavour indices, so does il(z). In the theory described 
by this action, the Euclidean space correlator of two baryon interpolating fields is given by: 

(0|Xa(a?)5Ca'(l/)|0>n = ^pj / VA^V$ Xa (x)x a '(y) ™P (-Sqcd^,^ - J d d z O(z) ■ Q(z)j (9) 



|n|->o 1 



Z[0] 



d d zn(z) f VA^V^ Xa (x)0(z)x a ,(y)e- s ^ +o(\n\ 2 ) 

- f d d z{l(z) (0\ Xa (x)0(z)x a ,(y)\0)o + 0(\n\ 2 ) 



\x a (x)x a '(y)\°)o 



Z[0] 

d d zn{z) (0|O(z)|0) 



where Z [Q] is the partition function of the modified theory, (■■■)□ indicates a correlator evaluated in this theory, and 
|fi| — > is a shorthand for requiring all tensor components of f2(z) to be small compared to Aq 1 ™^ 2 (from Eq. |(HJl, 

dimfi = d — dimO). For operators for which (0|O(z)|0) = [a notable exception being ip(z)ij)(z)], the (completely) 
disconnected contribution vanishes and we see that to leading order in the strength of the background field, O(z), the 
difference between the two-point correlators with fl(z) = and O(z) ^ is proportional to the form factors that we 
wish to extract. For example, defining the external field generalisation of Eq. (JSJ) as 

G% B (p', p, t;T) = J2 r 0a e- lp '- x e l P-y(O\xi B,) (x, t)xf\y, 0)|0) n , (10) 



then 



G BB (p,t;T) - G% B (p,p,t;T) t ^oo fi(0) tr [r {—ij> + M B ) A4bb [O] + M B ) 



Gbb (Pj t; r) |n|— o 2E B {p) tr[r(-^ + M B )\ 



(11) 



where fi(g) is the (discrete) Fourier transform of tl(z), M B is the ground state baryon mass, E B (p) = y/M B + \p\' 2 
and terms of 0(|f2| 2 ) are ignored. With appropriate choices of T, this gives the required form-factors at q 2 = 0. The 
non-forward case is discussed below. 

If O(z) depends on quark fields, the effects of the external field will manifest themselves in the quark determinant 1 
and valence quark propagators after the fermionic functional integration is performed in Eq. |j!JJ • In a lattice simulation 
where one approximates the integral over the field configurations by importance sampling, each choice of background 
field requires an additional ensemble of appropriately weighted gauge field configurations to be generated. We note that 
in quenched QCD a given set of gauge configurations can be modified to incorporate the effects of the exponentiated 
operator if it is composed of purely quark fields (such as the local electromagnetic and axial currents) as the operator 
and gauge field decouple in the absence of vacuum polarisation by dynamical sea quark loops. In more general cases, 
new gauge configurations are required even in the quenched theory. 

This exponentiated operator (external field) method is exactly what has been used to calculate the quenched 
magnetic moment of the proton in Refs. 0,0]. Here the addition to the Lagrangian is given by 

eB 



O em (z) = ^(zm(z), Sl^z) = — [7 , (12) 



1 One must be careful that the additional term in the action does not destroy the positivity of the determinant and thereby the probabilistic 
nature of the gauge integration measure. For example, adding the term ifi J d d z ip(z)~fQip{z) would be problematic. 
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producing a constant magnetic field in the Z3 direction over all lattice sites (ignore issues of periodicity at finite 
volume). In a lattice regularisation, the simplest local transcription of this current is not conserved and a multiplicative 
renormalisation factor is required; alternatively the lattice version of the conserved current [25L l27j can be used. For 
small field strengths, B, differences [T^.ll9|| between correlators of spin-up and spin-down baryons then determine the 
relevant magnetic moment, ^3, since the lowest energy eigenstates behave as Mg = Mb ± Hb\B\ + 0{\B\ 2 ) where 
± refers to spin anti-aligned or aligned with the magnetic field. To extract the magnetic moment, one requires one 
set of pure QCD gauge configurations generated with B = to determine the mass, Mb, and another independent 
set of configurations generated with the modified action with a\e B\ <C Aqcd to determine the shifted masses, Afg. 
In practice, a few ensembles of configurations generated with different values of B =/= may be needed to uniquely 
determine the mass shift linear in B. 

The quark angular momentum contribution to the proton spin can be measured from a forward matrix element [(| 
in a background field since it can be defined as a spatial moment, 3 q = J d 3 r r x T q (as with the magnetic moment, 
(j, = I f d 3 r r x j EM 28] ); in the background field approach 2 , one would add £1 times this integral to the QCD 
action and look at the shift in the exponential fall off of the two-point correlator exactly as for the magnetic moment. 
Since twist-two operators such as the energy-momentum tensor are not conserved even in the continuum, there is a 
multiplicative renormalisation factor that must be computed to obtain a final result. The renormalisation factors for 
various twist-two operators have been calculated |3(|. A more generic approach not relying on moment equations 
(and hence valid for more general operators that cannot be written in such a form) can also be used to determine J q . 
Nucleon matrix elements of the helicity independent, dimension four, twist-two quark operator take the form 

(N(p', s')\^ fl ^iD^ f \N(p, s)) = u N (p', s') [ 7 ^r } Af(q 2 ) - B f (q 2 ) + ^^C f (q 2 )] u N (p, s) , (13) 

where p = i (p + p 1 ) and {. . .} indicates symmetrisation of indices and subtraction of traces. In terms of the 
form- factors, Af, Bf and Cf, the total quark angular momentum content for flavour / is then given by 
Jf = \ [^/(0) + Bf(0)]. Since Bf(q 2 ) is always accompanied by a factor of the momentum transfer, the twist- 
two matrix elements must be calculated at q ^ (using Eq. J7J) for example) and then extrapolated to the forward 
limit. This necessarily leads to some uncertainty as the minimum available non-zero lattice momentum, q m i n — 27r/L, 
is set by the lattice size, L. For current lattice simulations, L < 3 fm, so q m i n > 0.4 GeV. Given such a large mo- 
mentum extrapolation, the use of the spatial moment definition in Eq. within the background field approach may 
provide the cleanest 3 determination of J q . 

Exponentiated operator methods can be used to calculate such off-forward matrix elements. To sec this, we consider 
the following correlator (p' ^ p) in the presence of an exponentiated operator J d d ztt(z) ■ O(z): 

Gg, fl (p / I p,t ; r) = ^r^e-p'-v^^oix^'^^Dxf^y^Dlo),, (14) 

= -EE r ^°( z ' T)e-P'- x e^(0| X (f'>(x, l)0(z, r)xf\y, -|)|0) + 0(\n\ 2 ) 

x,y z,r 

= - E E ^^T^r!(z,r) e -' i (P"-P')- x e l (P-P™)-y e -«(P m -P.)- Ze -(^+^)t/2 e (B^s m )r 

n,p n: s' m,p m ,s x,y z,r 

x (0\xi B 'H0)\n, p n , s')(n, p„, s'\O(0)\m, p m ,s)(m, Pm , S |^ S) (0)|0) + 0(\fl\ 2 ) 
= - E e HEn+EmW2 r 0a E e- lip - p ' yz e (E "- Em)T n(z, r) 

n,m s,s' z,r 

x (0\xi B,) (0)|n, P ', s') (n, p', S '|O(0)|m, p, s) (m, p, s\xf ] (0)|0) + 0{\n\ 2 ) 



t — >oo 



8 . 8 ' 



x( B ')(0)|i3'( P ^ S 0)(S'(p^ S 0|O(0)| J B( P , S ))< J B(p, S )|xf ) (0)|0), 



2 The difficulties of using continuum moment equations on the lattice highlighted in Ref. 29 ] do not apply to background field calculations. 

3 Another possibility is to consider the correlator subject to twisted boundary conditions 31,] (which correspond to a particular choice of 
external vector field) as they can reduce the minimum available momentum. 
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where fl(q) is again the (discrete) Fourier transform of fl(z), Eb and Eb> are the energies of the lowest B and B' states 
with momentum p and p'. Since energy and momentum are conserved in pure QCD simulations, G B , B {p', p, t; T) ~ 
5b'b3 3 (p' — p)- However, by including an external source that is inhomogeneous [fl = f2(z, t)] in the action, energy 
and/or momentum can be injected through the operator coupled to the source and off- forward correlators can have 
non-zero values. By choosing appropriate external fields, the form factors implicit in M.b'b[0] are easily determined 
from 



Gg, B (p',p,t;r') _ fi(?) tr[r'(-i/ + M Bl )M B 'B[0]{-ii> + M B )] 



Vg%b' (p'» p'> *; r) g bb (p, P , t, r) v*r [rFF + MbOI tr FH? + m b ) 



0m 2 ), (15) 



with relevant choices of T and T'. The simplest choice would be an external field that is a plane wave, ft(z) ~ 
exp(— iq ■ z), such that Sl(q) ~ <5 4 (<? — g), but other choices are possible. Correlators in which both the source and 
sink momentum are fixed are easily computed [3^ . 

This method will be applicable not only for off- forward matrix elements of the electromagnetic current (giving the 
Dirac and Pauli form factors) and the energy-momentum tensor (giving the form factors in Eq. 1|13[) ). but also for other 
quark bilinear operators. Important cases are the towers of twist-two operators O£ f" tln = ip jT^ L ° iD^ 1 . . .iD^ n 'ipf 
where = 7^, 7^75, cr Q/i and {. . .} again indicates symmetrisation of indices and subtraction of traces. The non- 
forward matrix elements of these operators correspond to moments of generalised parton distributions (see Ref. |3^| 
for a recent review). Finally we note that hadronic matrix elements of purely gluonic operators would also be 
accessible through the operator exponentiation approach but operators involving more than two quark fields cannot 
be exponentiated as they would prevent the fermionic functional integrals from being integrated exactly. 

For many quantities the background field approach would not be particularly appealing as the added overhead (with 
respect to calculating a baryon two point correlator) of computing the quark propagators required in the standard 
three point correlator analysis is minimal as compared to the cost of generating additional ensembles of dynamical 
gauge configurations; particular examples are the moments of the isovector parton distribution functions. However, 
for flavour singlet quantities the situation may be reversed. In the operator insertion approach, one must evaluate 
the contributions of quark-line disconnected diagrams which are very difficult to compute as they involve propagators 
from all points on the lattice to themselves. Consequently, they require ensembles of gauge configurations that are an 
order of magnitude larger to achieve the same precision as in the corresponding connected diagrams [T^ | even when 
stochastic estimator techniques 0, ES HE E3> EE (which improve greatly on the brute force approach) are applied. 
In the background field approach, disconnected contributions are included automatically; for these calculations one 
requires a few additional, moderately sized ensembles of configurations generated with modified actions and this 
may be a more efficient means of calculating flavour singlet quantities. In terms of statistics and identification of 
the ground state, the external field method is somewhat similar to the plateau accumulation method [3!|, but with 
the added advantage of automatically including quark-line disconnected diagrams. To extract physical quantities in 
either approach, calculations must be repeated at different quark masses, volumes and lattice spacings and then the 
necessary extrapolations must be performed. This may be more difficult using external fields than in the operator 
insertion approach where multi-mass techniques |4£j are well developed (though similar techniques may also make the 
background field calculations easier). Numerical work to investigate the relative costs of the different approaches is 
encouraged. 

Two baryons: In the two baryon sector, exponentiation of operators will also prove useful. It has recently been 
proposed that background electroweak fields can be used to measure the electroweak properties of two-nucleon states 
such as the deuteron magnetic moments and the cross-section for neutrino breakup of the deuteron p^ |. For such 
two-hadron systems, even a calculation of the lowest energy levels involves large numbers of distinct Wick contractions 
and is only presently becoming feasible |4l| ; calculating matrix elements of such systems by an operator insertion is 
probably beyond the limits of current computational power. However, by measuring finite volume shifts in the low-lying 
two-particle energy levels in appropriate background fields, the various electroweak properties can be determined using 
effective field theory. Again, the exponentiated operator approach is especially suitable for flavour singlet quantities 
as disconnected diagrams, which would be required in the operator insertion approach, are included automatically. 
Provided one can construct deuteron interpolating fields that sufficiently project onto its different polarisations, similar 
background field methods will yield the quadrupole moment of the deuteron. 

Another important phenomenon which background fields may prove useful in investigating is the modification of 
structure functions in nuclei - the EMC effect [16J . The difference of the ratio F$(x)/ [F$(x) + F£{x)\ from unity is the 
simplest quantity in this regard (though experimentally the magnitude of the effect increases with atomic number). 
Important matrix elements that one might consider to investigate this on the lattice are {N N\O^ ^^\N N) [l5j 
since these determine the structure functions of two-nucleon states such as the deuteron via the operator product 
expansion. In a similar manner to the deuteron magnetic moment and neutrino-deuteron breakup cross-sections |24| , 
one can exponentiate the twist-two operators to create background fields. Then evaluating two-particle energy levels 
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in these field and matching to effective field theory will enable calculation of the required matrix element. Again, 
since these matrix elements are isoscalar, the external field approach should be competitive with operator insertion 
techniques 4 . Finally, we note that the same sets of generalised background field gauge configurations that one might 
use to calculate moments of twist-two operators in the single baryon and meson sectors would also be appropriate 
here. 

In summary, background field methods in which an operator is exponentiated into the action are powerful tools for 
evaluating general hadronic matrix elements in lattice QCD. Here, we have highlighted some of the features of this 
approach. The method is not restricted to operators whose background fields have classical analogues; for example, 
twist-two operators, such as that determining the total angular momentum content of the proton, are ideal candidates. 
For flavour singlet matrix elements, disconnected diagrams are automatically calculated in the external field procedure, 
potentially leading to increased statistical precision in these quantities. The technique is also not limited to forward 
matrix elements; objects such as the electromagnetic form-factors and generalised parton distributions can also be 
computed with background fields. These advantages come at the cost of requiring additional ensembles of gauge 
configurations to be generated. 

The author is grateful for discussions with C.-J. D. Lin, W. Mclnitchouk, M. J. Savage, A. Shindlcr and J. M. Zanotti. 
This work is supported by the US Department of Energy under contract DE-FG03-97ER41014. 
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